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Abstract 

We discuss the various aspects of two-dimensional QCD2{N — > oo) 
(the 't Hooft model[[l|]). Our main interest (motivated by the corre- 
sponding analysis in the four dimensional QCD) is the vacuum struc- 
ture of the theory. We use the very general methods in the analysis, 
such as dispersion relations and duality in order to relate the known 
spectrum of QCD2 to the different vacuum characteristics. 

We explicitly calculate (in terms of physical parameters like masses 
and matrix elements) the chiral condensate as well as the mixed vac- 
uum condensates: 

(0\q(9e^G^) n q\0) ~ M% f (0\qq\0). 

We interpret the factorization property for the mixed vacuum con- 
densates as a reminiscent of the master field at large N. 
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1. Introduction 



In this paper we analyze the vacuum properties of two-dimensional QCD. 
The standard way of doing so is the "canonical" (or non-canonical) gauge 
fixing, analysis of different constraints, elimination of redundant degrees of 
freedom e t.c. We take the opposite approach. We assume that the spectrum 
is known and our main goal is find out the vacuum properties which would 
correspond to the given spectrum. 

Since the pioneering paper by 't Hooft|T| in 1974, QCD 2 (N — > oo) has 
been subject of many investigations @ - 0- The list of problems discussed 
in these papers is impressive: scaling and e + e~ annihilation; deep-inelastic 
scattering and Drell-Yan formula; Regge behavior, analysis of form factor 
and wide angle behavior of exclusive amplitudes; light cone quantization and 
e t.c. 

The next step which has been taken is the analysis of the vacuum prop- 
erties which would correspond to the well established spectrum. It has been 
realized that the vacuum of the theory in the 't Hooft limit 

q 2 N ~ const. N — > oo, m a ^> q ~ — = (1) 

q JN 



g 2 N 
12n 



is quite nontrivial. In particular, the chiral condensate (qq) = —N 
in the limit m q —> has been calculated ||0.The result was confirmed by 
numerical and independent analytical calculations [IT|. Moreover, the 

method has been generalized for the nonzero quark mass and the correspond- 
ing explicit formula for the chiral condensate (qq) with arbitrary m q has been 
obtained fL2| . Let us note that there is no contradiction with the Coleman 



theorem |13| at this point, because the BKT (Berezinski-Kosterlitz-Thouless) 
behavior takes place in the large N limit [Q. 

Recently there has been a renewal of interest in the study of QCD 2 Hl4| - 



The motivation for this interest was quite different for different people: 
it ranged from an attempt to find a string representation of QCD in four 
dimensions to mastering the instantons and master field in this theory. The 
important lesson to be learned from the new development, we believe, can 
be formulated as follows: the structure of gauge theories even in two dimen- 
sions (with and without matter fields) is very complicated. However, some 
physical questions can be formulated irrespective of number of dimensions. 
In certain circumstances , the correctly formulated questions can even be 
answered. This gives some hope (and hints) that the similar formulation of 
the analogous problems hopefully can be found in four dimensions. 

This paper is largely motivated by analysis of nonperturbative wave func- 
tions with a minimal number of constituents in four dimensional QCD. As 
is known, such a function gives parametrically leading contributions to hard 



3 The chiral limit m q — > can be considered only after the limit N — > oo is taken. These 
limits do not commute, see the next section for explanation. 
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exclusive processes. The corresponding wave functions within QCD have 
been introduced into the theory in the late seventies and early eighties [33 
to describe the exclusive processes. We refer to the review papers |35|[, [j36[ 



on this subject for the detail definitions and discussions in the given context 
of Wilson operator expansion. 

The two-dimensional QCD 2 (N — > oo) is the ideal toy model for such 
kind of problems, because the only physical states in the model are states with 
minimal number of constituents (the qq mesons)^. The properties of these, 
nonperturbative wave functions might be quite nontrivial and unexpected. 
It would be very interesting to check the analytical nonperturbative methods 
(mainly the QCD sum rules) which have been developed for such an analysis. 

As is known, the QCD sum rules approach operates with objects like vac- 
uum condensates of different local operators. Eventually, if we knew all types 
of vacuum condensates, we would calculate an arbitrary correlation function 
and thus, through the dispersion relation, we would find the spectrum and 
amplitudes. 

The subject of the present paper is analysis of the vacuum structure. 
More specific, we are interested in the exact calculation of the chiral, mixed 
vacuum condensates (0|g((?e Ai jyG ! ^) n g|0). Once these nonperturbative con- 
densates are found, than some more complicated nonlocal vacuum expecta- 
tion values (like Wilson line (0|W(a:)|0) = {0\q{x)e ig fo A » dx »q(0)\0) and its 
superpositions {0\W{xi)W{x 2 )W{x 3 )...\0)) can be in principle evaluated. 

The paper is organized as follows. In Sect. 2 we review some QCD 2 (N — > 
oo) properties with emphasis on the vacuum structure of the theory. We 
review the spectrum found by 't Hooft and explain why this spectrum un- 
avoidably leads to the existence of the chiral condensate {qq). Through the 
dispersion relation we explicitly express the condensate in terms of spectrum 
and physical matrix elements. We qualitatively explain why the nonzero 
magnitude for {qq) does not contradict to the Coleman theorem and why the 
chiral limit m q — > can be considered only after the large N limit is taken. In 
different words, we demonstrate that to preserve the 't Hooft solution (which 
corresponds to the selection only planar, leading at large N, diagrams) we 
need to require that the inequality m q ^> g ~ l/y/N is fulfilled. Otherwise, 
the nonplanar diagrams which might have the singular factors like 1 / m q come 
into the game and destroy the 't Hooft solution. 

Sect. 3 is devoted to the calculation of the mixed vacuum condensates. 
We try to avoid, in all cases, the discussion of such complicated problems as 
infrared and ultraviolet regularization, the problem of splitting of two local 
operators, and many other problems which unavoidably accompany any cal- 
culation in local field theory. Instead, we use the knowledge of existence of the 
nonperturbative vacuum condensate {qq) to calculate the more complicated, 
higher dimensional condensates^. The crucial reason for such incredible 

4 Let us remind that an additional creation of qq pairs is suppressed by a factor 1/N. 

5 The gluon in QCD2 is not a physical degree of freedom. Thus, one could anticipate 
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simplification is the factorization property (Qi ■ Q 2 ) = (Qi)(Q2) + 0(1/N) of 
vacuum expectation values in the large N limit. Thus, all vacuum conden- 
sates can be reduced to (qq) and its powers. We believe that if the master 
field is found, it should reproduce the vacuum properties mentioned above. 

Sect. 4 is our conclusion and outlook. To be more specific, we discuss 
some phenomenological consequences of the obtained formulae which might 
be interesting for the analysis of nonperturbative wave functions and heavy- 
light quark system. 

2. The spectrum and chiral condensate in QCD 2 (N — > oo). 

The model we shall consider consists of quark in fundamental represen- 
tation interacting via an SU(N) color gauge group. We follow the notation 
of ref.|2| and present the 't Hooft equation in the following form: 

where symbol P notes as the principal value of the integral, and < x < 1 
is the fraction of the total momentum of the bound state carried by quark q 
with mass m q . The quantity m\ = 3 -^- is the basic mass scale in the theory 
and the index n classifies the ordering number of the bound states \n,p) with 
total momentum p^. The same wave function can be expressed in terms of 
the following matrix element ||: 

M%) =f^f dy + e- iy+{1 - 2x)p -(0\q(-y)q(y)\n,p)\ y ^ =0 . (3) 

Let us note that matrix element on the right hand side is written in the light 
cone gauge A_ = 0; to restore the manifest gauge invar iance one can insert 
the standard exponential factor e tg -I ' A ~ dy + into the formula (0). 

Let us review some important properties of equation (0). The entire spec- 
trum is discrete and classified by the integer number n. The wave functions 
4> n {x) are orthogonal, complete and obey the following boundary conditions 

9 

(j) n (x) -»• [ar(l - x)f, x^O, x -»• 1, it (3 cot (it (3) = 1 1. (4) 

It turns out the spectrum of states is almost linear for large n: 

m 2 n ~ TT 2 mln, (f) n (x) ~ v^2sin(7rnx) (5) 

and does not depend on mass of the quark. What is more important, in the 
chiral limit m q — > the lowest level (we call it 7r meson) tends to zero m 2 ~ 
77i q and one could expect the nonzero magnitude for the chiral condensate. 

that any local operator can be expressed in terms of the quark fields. As is expected, this 
is indeed the case. 
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We define the chiral condensate in the current algebra terms as follows: 



= lim i / d 2 xe tpx d^\T{q ltllh q{x), q^W) = (6) 

PfJ. *V J 

2z(0|gg|0) + 2m q ■ (0|T{gz 75 g(x), gz7 5 g(0)}|0). 

As we already mentioned, the only states of 't Hooft's solution are the quark- 
antiquark bound states. Thus, they must saturate the dispersion relation. 
Upon inserting this complete set of mesons to the (§) one thus obtains: 

<0|g#> = -m g £^-§, (7) 
where f n is defined in terms of the following matrix elements 

/nl- I \ t f m i f 1 a (a\ 

(0| TO |n) = S^-fn, fn=^-J Q Jy^** (8) 

In the chiral limit the only state which can contribute to the formula (0) is 
the 7r meson. Its matrix element can be calculated exactly and we end up 
with the following expression for the chiral condensate in the m q — > limit 



(0|gg|0) = —N-^=, ml = f n=0 = m\ = m q ■ (9) 



As was expected, we find that (0|gg|0) ~ N. Besides that, as we already 
noticed in |J, if we put m q = from the very beginning, then (0|gg|0) = 0. 
This corresponds to the different regime when m q <C g ~ \/\fN, when 
nonplanar diagrams come into the game. We discuss this point a little bit 
later. The last remark is the observation that the entire nonzero answer for 
the condensate comes from the infrared region of the integration in eq.([8]): 
x ~ 0, x ~ 1 which corresponds to the situation when one of the quarks 
carries all the momentum and the second one is at rest. 

The sum (0) can be calculated exactly for arbitrary m q |T^| . The crucial 
point is that for arbitrary m q the nonzero contribution comes from the highly 
excited states {n ^> 1) only. The properties of these states are well-known: 

j?2 2 2 2 2 2 1 /in\ 

fn -» 7T m , m n -> 7T m ■ n, n > 1, (10) 
and thus the sum (|7|) can be explicitly evaluated with the result |L2] : 



m N 1 
(0|gg|0> = -^-{log(Tra) - 1 - lE + (1 - -)[/(«) - al(a) - log 4]}, (11) 

III a 



where a = 7^ = 0.5772.. is Euler's constant and 

dy 1 

sinh y cosh y 



y 



1(a) 



y 2 [a(y cothy — 1) + 1] 
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This result is exact for large N and arbitrary quark mass as far as we remain 
in the 't Hooft regime ([[]), i.e. m q ^> g ~ 1/yN. It reduces to the eq.(Q) in 
the limit a — > oo, as it should. 

The last condition (m q ^> g ) which has to be satisfied for the 't Hooft 
solution to be valid, requires some additional explanation. Roughly speaking, 
nonplanar diagrams may contain a factor ~ m~ x which at m q = blow up 
and the theory changes completely. The concept of the proof that there exists 
a factor ~ m" 1 in nonplanar diagrams is the following. 

Let us consider the correlation function for p — > 

i J d 2 xe tpx (0\T{qq(x), gg(0)}|0) = P(p 2 ) (12) 

The 't Hooft solution suggests that only planar graphs are taken into account 
and, consequently, the spectral density contains only the contribution of one 
meson states. For these contributions P p i an ar ~ N. In the chiral limit, we 
can calculate the two-pion contribution exactly! This contribution is not 
accounted for in deriving (0). Of course, the two-pion contribution is 

2 

suppressed by a factor 1/N. However, it contains a term ~ which tends 
to infinity for m q — > 0. The presence of the factor ~ m" 1 in nonplanar 
diagrams leads to the aforementioned constraint on m q (|I|). 

Now, let us explicitly demonstrate the existence of the term ~ m" 1 for 
the two-pion contribution. In order to do so, let us write down a dispersion 
relation for P: 

P(0) = - / —ImP{s), (13) 

7T J Ami s 

where ImP(s) is the physical spectral density. The -kit contribution is fixed 
uniquely by (|9|), because of the special role of pions ||: 

tttt qq Uo = -ImP^(s) = -± 14 

V3 vr 6 J s{s -Ami) 

P™(0) = 



mlir 2 r°° ds m^Ti 2 1 



6 J^l S J s(s -Ami) 12m l m g 



It is clear that the only cause for a singular ~ 1/ m q behavior is the finiteness 
of the pion matrix elements at zero momentum. At the same time this 
contribution does not contain the large factor N which accompany a one 
meson contribution to the same correlator. To suppress these nonplanar 
diagrams we need to require iV 3> ^f. Thus, we would expect that some 
kind of phase transition may occur in the region m q ~ g, where we would 
expect a complete restructuring of the theory. 

The last subject we would like to discuss in this section is the strict 



Coleman theorem [13| which states that a continuous symmetry cannot be 
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broken spontaneously in two dimensional theories. As we discussed earlier 
H] we expect that as in the SU(N — > oo) Thirring model (where the chiral 



symmetry is "almost" spontaneously broken f37|), the BKT effect |3§] op- 
erates in regime (|ip. This fact also confirms the 't Hooft spectrum: states 
with opposite P parity are not degenerate in mass and there is an "almost" 
Goldstone boson with ~ m q + 1/N. 

To be more specific, one can show || that in QCD2(N — > oo) the behavior 
of the proper two-point correlation function is as follows: 

(0\T{q L q R (x), q R q L (0)}\0) ~ x~k . (15) 

Such a behavior together with claster property at x — > oo implies the ex- 
istence of the condensate at N = oo in a full agreement with our previous 
discussion. At the same time, for any finite but large N, the correlator falls 
off very slowly demonstrating the BKT-behavior with no signs of contradic- 
tion to the Coleman theorem. 

3. The vacuum expectation value of the mixed vacuum conden- 
sates. 

We start from the consideration of the simplest mixed vacuum expectation 
value 

(qD^D u q) = ^^(qDxD^xl* ~ eAa7s)<?) = ^^(qmxvGxalsq) , (16) 

where [_D M , D v ] = —igG^ is the field strength tensor of the gauge field and 
we have used the well-known properties of the 7 M matrices in two dimensions. 
Thus, in the chiral limit, m q — > 0, the operator we are interested in can be 
expressed exclusively in terms of the field strength tensor G^ v . 

What's more, in two dimensions an arbitrary operator can be reduced to 
the form which contains the field strength tensor ige^G^ only. Indeed, the 
covariant derivatives (qD n ■ D^D^q) placed at the very right and at the very 
left (near the quark fields) can be transformed into the operator ige^G^ as 
before. To do the same thing with operators which placed somewhere in 
the middle, we need to act (say) on the right until the quark field is reached. 
By doing so, step by step, we create many additional terms which are either: 
commutator like [_D M , D u ] = —igG^ which is the field strength operator or 
commutators like ~ [D\, e^G^]. Fortunately, in two dimensions these terms 
are related to creation of the quark- antiquark fields and we discard them in 
according to large iV counting rule 

Indeed, the quark-antiquark operator comes from the equation of motion 
DpG^y ~ q^ u q when we sum up over the common index \i. This is true in any 
number of dimensions. In two dimensional case, even without a common in- 
dex fx, the action of on field strength tensor produces the quark-antiquark 

6 Of course this is not the case in four dimensions where [D^Gxa-] is independent 
operator which can not be reduced to some quark fields. 
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operator. This can be seen by noting that in 2 dimensions G^ v {x) ~ e lu/ E(x) 
with scalar function E ~ cao-Gao- and D^E ~ qj^Jsq- 

Thus, we end up with the vacuum condensates (q(gE) n q) which are ex- 
pressed exclusively in terms of the field strength tensor and our nearest prob- 
lem is their calculation. 

Here we sketch the idea of this calculation. We choose the light cone 
gauge, v4_ = -j^(A — Ai) = 0. In this gauge we have the usual constraint 



in the gauge sector (Gauss law) (for more details see e.g. [Ill 



d_E ab ~g( q V> q l-±-5 ab ql c ql), 

where a, b are the color indices. Here the right moving fermion q + are dy- 
namical degrees of freedom; the left-moving fermion q_ are non-dynamical 
degrees of freedom in this gauge. The latter can be eliminated by the follow- 
ing constraint: 

<9_g_ ~ m q q + . 

The next step is to use the Gauss law to calculate the mixed vacuum 
condensate. For the simplest case it looks as follows: 

(qige^G^q) ~ (ql a gE ab q b + ) ~ (ql a g 2 i(£g?)g*> 

Although the explicit expression of the Green function -g- is well known (it 

is the step function e(x_ — x'_) in the coordinate space, see e.g.|ll[), as will 
be clear soon, we do not need its manifest form. 

What does matter, is the important property of the large N limit that 
reduces the expectation value of a product of any invariant operators to 
their factorized values Thus, the condensate under consideration can be 
presented in the following form 

(qige»„G^ 5 q(0)} ~ ^(O)^)) ■ (^-gj?£(0)) 

The first term on the right hand side of this expression is reduced to the chiral 

X 71 

condensate (qq). Indeed, the additional terms like -^{q{d~) n q) that come 
from the Taylor expansion, are zero because of the Lorentz invariance of the 
vacuum state. The second term is also reduced to the chiral condensate if 
one takes into account the constraint mentioned above and which we would 
like to write down in the following way: 4-q+ ~ ^- q~ 



As we have already mentioned, neither explicit formula for the Green 
function has been used, nor specific regularization has been assumed in 
this derivation. 

The final formula for the simplest mixed condensate (0) which accounts 
for all numerical factors takes the following form: 

= (17) 
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Few comments are in order. First of all, as we expected, the higher 
dimensional condensate can be expressed in terms of the fundamental chiral 
condensate (plfFID . We note also that the mixed condensate (|17D has the same 
N dependence as the fundamental one, (qq) (remember, g 2 (qq) ~ g 2 N ~ 1). 

The sign of the right hand side in eq.(|TT|) is the result of the calculation. 
We would like to pause here to compare this result with somewhat analogous 
in four- dimensional QCD, where the corresponding formula looks as follows: 

±(q(x,Dj 2 q) = lx 2 • (qiga ^G^q) ~ ± x 2 0.8GeV 2 (qq) . (18) 
In both cases, the sign of the ratio 

(q(x,D,) 2 q) x2 (qDjq) ^2 < Q 
(qq) (qq) 

is negative for negative x 2 , where operator expansion is effective. In four 
dimensions this sign has very deep physical meaning, because there is one- 
to-one correspondence between such kind of the vacuum condensates and 
mean values of transverse moments for the pion wave function POH 7 !- There is 
no such interpretation in 2d (no transverse direction in this case), however, 
one can argue that for the space- like interval, the sign in the eq.(|17D does 
correspond to the positivity of the Hermitian operator (iDi) 2 . 

The next step is the application of the same procedure for the operators 
with arbitrary number of gluon field insertion E ab . It can be done in the same 
way as before just because of the factorizability mentioned above. Thus we 
arrive to the following formula 



(q{DM n q) = ^(q{ige x<J G Xal sYq) = (-£M)»<w), (19) 



which is our main result. The comment to this formula can be formulated in 
the following way. Each time the insertion of an additional factor (D^D^), 
which is proportional to the field strength tensor gE, gives one and the 
same numerical factor (|19|). Situation can be interpreted as we would have 
a classical master field |39j which we insert in place of gE in the vacuum 
condensates. Because of its classicality, it gives one and the same numerical 
factor. 

The second important msg. is as follows. The vacuum condensate of 
an arbitrary local operator can be reduced through the equation of motion 
and constraints to the fundamental quark condensate (§,|11]). This is exactly 
what one could expect in QCD 2 where a gluon is not a physical degree of 
freedom, but rather is constrained auxiliary field which can and should be 
expressed in terms of quark fields. 



7 This is not a big surprise, however, because the pion is Goldstone particle and its 
matrix elements very often can be reduced to the vacuum condensates through the PCAC. 
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With these remarks in mind one could calculate (in principle) the vac- 
uum expectation value of an arbitarry nonlocal operator like string operator 
(W(x)) and its superpositions (W(xi),W(x2),W(x3)...). The only what we 
need to do is to use the Taylor expansion 

n=oo i 

(WO = (0|g(*)Pe^o A ^q(0)\0) = £ Mx,D,f n q) (20) 

n=0 \ ln >- 

for the nonlocal operator we are interested in and substitute the values for 
the vacuum condensates we calculated previously ([19D . However we better 
stop here. 

5. Conclusion and Outlook. 

We have calculated the mixed vacuum condensates (|I^) and have ex- 
plained of how to reduce an arbitrary local vacuum condensate to already 
known values (|19|) . We also have demonstrated of how to calculate the vac- 
uum expectation value of the nonlocal operators ( like string operator and 
its superpositions) by mean the Taylor expansion and reducing the problem 
to the previous one flUf). We shall not discuss here the theoretical issue of 
related problems (which are quite interesting, by the way). Instead, in this 
conclusion we would like to mention some "phenomenological" applications 
which might be interesting for the future investigation. 

First of all, as we already noticed, (W) naturally appears in the analysis of 
the heavy-light quark system qQ within operator product expansion. Indeed, 
if we consider along with f4lJ,J|2| the correlation function (T{qQ(x), Qq(0)}), 



describing this system, we end up (in the limit Mq — > oo) with the object 
which completely factorized (in accordance with HQET^3|) from the heavy 
quark and which was called NLC (nonlocal condensate) 1 



(T{qQ(x),Qq(0)}} ~ (q(x)P exp(ig / A^dx^q^S)) + perturb, part. (21) 



o 



All nontrivial, large distance physics of the system is hidden there. It has 
been noticed recently that this system might provide a definition of the con- 
stituent quark in QCD. Together with perturbative contribution one should 



expect the following behavior for this correlator] 42 1: 



(T{qQ(x),Qq(0)})~ e - A -*. (22) 

We do not expect that the function (W) can provide such a behavior by itself 
in t'Hooft model. The reason for that is related to the perturbative terms, 
proportional to ~ ( a_ )" (x 2 ) n . These contributions go on the same foot as 
nonperturbative ones due to the dimensional coupling constant g in 2d, and 
they interfere with expansion (|20|). 

Naively, one could expect that nothing like that might happen in four 
dimensional QCD4, where the coupling constant is dimensionless and no 
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power corrections might occur in perturbative series. However this is not 
completely true, because of the so-called renormalons, which may provide 



some effective power corrections ||44j| . 



Another, but related issue is as follows. As is known, the analysis of 
different problems (like inclusive decays, distribution function...) in heavy- 
light quark system, requires the precise information about the corresponding 
behavior in the so-called end-point domain. Formally, it requires the knowl- 
edge of the infinite set of matrix elements like (Hq\QD ^....D n n Q\Ho) . Our 
remark is that such matrix elements can be explicitly calculated in QCD 2 
as we discussed before. We believe that QCD 2 as a toy model may provide 
at least a hint on analytical properties of this, so-called universal distribu- 
tion function F(y) This universal distribution function can be presented as 
matrix element F(y) ~ / e ixy (H Q \Q(x)e ig So A » dx »Q(0)\H Q ). Operator ex- 
pansion presents this (presumably regular function F(y)) as a series of delta 
function and its derivatives F(y) ~ S«n^"(l — y)- Thus, to find out the 
analytical properties of such a function in the real world is a difficult task 
which requires the knowledge of distant terms of the series. The lessons 
which might be provided by QCD 2 could be useful. 

At the end, then, we come back to the beginning. We have already 
noticed in the Introduction that QCD 2 is the nice model for analysis of 
nonperturbative wave functions themselves and methods which one can use 
to extract the corresponding information. We refer the reader to the recent 
papers on this subject [46|,|47] for the references and new development. Here 



we would like to note that in QCD 2 we know the nonperturbative wave 
functions as well as whole set of local condensates (0) which play the crucial 
role in this study. Thus, one may try to check the methods which have been 
developed for corresponding analysis. 

It is a pleasure to thank Professors Gary McCartor and Kent Hornbostel 
for encouragement and interesting discussions related to the subject of this 
paper. 
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